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ABSTRACT 

We present interferometric observations of the Be star ( Tau obtained using 
the MIRC beam combiner at the CHARA Array. We resolved the disk during 
four epochs in 2007-2009. We fit the data with a geometric model to characterize 
the circumstellar disk as a skewed elliptical Gaussian and the central Be star as 
a uniform disk. The visibilities reveal a nearly edge-on disk with a FWHM 
major axis of ~ 1.8 mas in the if-band. The non-zero closure phases indicate 
an asymmetry within the disk. Interestingly, when combining our results with 
previously published interferometric observations of ( Tau, we find a correlation 
between the position angle of the disk and the spectroscopic V/ R ratio, suggesting 
that the tilt of the disk is precessing. This work is part of a multi-year monitoring 
campaign to investigate the development and outward motion of asymmetric 
structures in the disks of Be stars. 

Subject headings: circumstellar matter — stars: emission-line, Be — stars: indi- 
vidual {( Tau) — techniques: interferometric 
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Introduction 



C Tau (HR 1910, HD 37202, HIP 26451) is a bright Be star with an extensive observa- 
tional history that spans photometric, spectroscopic, polarimetric, and interferometric tech- 
niques. The spectrum of C Tau shows the standard dou ble-peaked Ha profile, indicative of a 
disk in Keplerian rotation (e.g. iPorter &: Riviniusll2003[ ). The relative height of the blue- and 



red-shifted peaks of t he emission lines {V/R ratio) shows cyclic variab i lity measured on the 
order of ~ 1429 days jRivinius et al.lbooel : IPoUmann k RiviniusI 120081 iRuzdiak etaP 12009 



Stefl et al.l 120091 ). The spectra of C Tau also occasionally show more complex triple-peaked 
and "shell" profiles. Current models suggest that the line profile variations in C Tau and 



other Be stars can be explained by global one-armed oscillation models (e.g. ICarciofi et al. 
2009h . 



C Tau is also a single- lined spectroscopic binary with a period of 133 days (e.g. iHarmanec 



1984J : iRuzdjak et al.ll2009l . and references therein). The companion has not yet been detected 
directly. By placing limits on the orbital inclination, the mas s function indicates that the 
companion is a low-mass object(~ 1 M©; iRuzdjak et al.ll2009l ). several magnitudes fainter 



than the primary Be star. The companion could be a main seq uence star, a neu tron star, a 
w hite dwarf, or an evolved hot subdwarf like that found in Per (IGies et al.lll998l ). However, 



as 



Ruzdjak et al.l (120091 ) discuss, a hot subdwarf may heat the outer facing regions of the 



disk, producing narrow emission lines like He I A6678. The absence of such e mission in the 



spectr um of C Tau suggests that the companion is not a strong fiux source. iFloquet et al. 



( 119891 ) rule out a cool luminous giant based on measurements of the infrared fiux. 



Given its brightness (Knean=3.0, ifmean=3.0) and distance ((i=128 pc; iPerryman et al. 



19971 ). C Tau is an ideal source for studying the structure and dynamics of its circumstellar 
disk using optical/infrared (IR) interferom etry. Early spatially resolved measurements re- 
vealed the elliptical shape of the C, Tau disk ( iQuirrenbach et al.lll994j . 119971 : [Baldwin fc Hanifi 
2OO2I : iTycner et al.ll2004j ). indicating that the circumstel lar materia l is con tained in a fiat- 
tened disk inclined nearly edge-on to the line of sight. iGies et al.l (12007! ) determined the 
geometry and density structure by fitting an isothermal model of a disk in Keplerian ro- 
tation to CHARA Classic interferometric data of the i^'-band emission. Using the GI2T 



interferometer, IVakili et al.l (119981 ) detected an asymmetry in the disk of C, Tau based on 
differential phases measured across the Ha line. These results suggest a shift in the position 
of a bulge located in the disk, consistent with the prograde motion of a one-armed spiral 
oscillation in the disk. Moreover, the nature of the asy mmetry and its rel ation to the tempo 



ral V/ R profile variations were investigated further by IStefi et al.l (120091 ) and ICarciofi et al. 



( I2OO9I ) using VLTI/AMBER observations. The visibilities and differential phases measured 
across the Br7 line are consistent with an oscillation pattern created by a one-armed spiral 
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in the disk. 

In this paper we present muhi-epoch interferometric observations of ( Tau obtained 
using the CHARA Array in the if-band. The visibihties provide information about the size 
and orientation of the disk while the closure phases indicate the presence of an asymmetry 
within the light distribution. We fit the data from each epoch using a model consisting of 
the central star and a skewed elliptical Gaussian disk. We present the results in § [31 In §U 
we describe the physical characteristics of the model and discuss the changes we observe in 
the orientation of the disk on the sky, the asymmetric light distribution, and their relation 
to the cyclic variations measured in the emission line profiles. In §4.2[ we outline the details 
of a precession model that can explain the changes we observe in the disk of ( Tau. We 
summarize the results of the study in §3 



2. CHARA Array Observations 



The CHARA Array is an optical/IR interferometer located on Mount Wilson (jten Brummelaar et al. 



20051 ) . The array has six 1-meter telescopes arranged in a Y-configuration with baselines 



ranging from 34-331 meters. We used the Michigan Infrared Combiner (MIRC: lMonnier et al 



2004J . l2006af ) at the CHARA Array to observe the disk of C Tau in the if-band. MIRC com- 
bines the light from four telescopes simultaneously, providing visibility amplitudes on six 
baselines and closure phases on four triangles. It uses single-mode optical fibers to spatially 
filter the light. The fibers are brought together by a V- groove array in a nonredundant 
pattern which encodes the overlapping fringes formed from the outgoing light with distinct, 
spatial interference frequencies. We used the low spectral resolution prism (R~50) to dis- 
perse the fringes across eight spectral channels in the H-hand (A = 1.5-1.8 fim). 

Table [1] provides a log of the observations that lists the UT date of the observation, the 
configuration of telescopes used, and the observed calibrator stars. On 2007 Nov 11-14 we 
used the inner array (S2-E2-W1-W2) of shorter baselines ranging from 108-248 m. On the re- 
maining dates we used the outer array (S1-E1-W1-W2) of longer baselines ranging from 108- 



lar diameters smaller than 0.9 mas (selected from the catalogs of 


Pasinetti-Fracassini et al. 


2001 




Merand et al. 


2005 




Richichi et al. 


2005 




van Belle et al. 


2008 


). We derived angular 



diameters of the calibrators by fitting all the available flux data with reddened models of 
the spectral energy distribution. In the best cases, we included UV-fiuxes from the Inter- 
national Ultraviolet Explorer satellite, optical spectrophotometry, and near-IR fluxes based 
upon 2MASS and other sources. The more poorly studied stars have flux estimates at least 
for the Johnson UBV bands and for the JHKg 2MASS bands. The final ingredient in the 
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fit is an estimate of the stellar effective temperature Tes that usually comes from detailed 
spectroscopic investigations. The observed spectral energy distributions are then fit with 
low resolution spectra from the models of R. Kurucjil (for T^s > 5000 K) or Gustafsson et 
al. (2008; for Teg < 5000 K). Table |5] lists the names and HD numbers of the calibrators, 
spectral classification, V and H magnitudes, the adopted Tes and reference source (also for 
the adopted value of gravity logg), as well as the derived interstell ar reddening (usually for 
a ratio of total to selective extinction, R = 3.1; iFitzpatrickl Il999l ) . and the limb darkened 
angular diameter. In what follows we assume for simplicity that the limb darkened angular 
diameter equals the uniform disk angular diameter (since the differences are very small in 
the near-IR where limb darkening is minimal). 



The data were reduced using the standard MIRC reduction pipeline (e.g. iMonnier et al. 



20071 ). The reduction process involves Fourier transforming the raw background-subtracted 
data to obtain fringe amplitudes and phases. Corrections are then applied to the visibility 
amplitudes to account for the fiber coupling efficiencies to correct for the different flux levels 
in each of the four beams. The data from 2009 Nov 10 were acquired using the photometric 
channels recently installed in MIRC that divert 20% of the li ght going through each fiber to 



measure directly the contribution of light from each telescope (IMonnier et al.ll2008l : IChe et al. 



20101 ). Drifts in the overall system response are calibrated by using single stars of known sizes 
observed before and after the target observations. The reducti on pipehne outpu ts calibrated 
squared visibilities and closure phases in the OIFITS format@ (jPauls et al.ll2005l ). The errors 
in these values are derived by combining the scatter in the measured data with the errors 
propagated through the calibration process. 

Figure [1] shows the u-v coverage on the sky sampled by the CHARA Array during the 
epochs of the ( Tau observations obtained with MIRC. The squared visibilities measured on 
the six baselines are shown in Figure HI The range of angular sizes indicated by the visibilites 
along different projection angles reflects the elliptical shape of the ( Tau disk on the plane of 
the sky. The closure phases on the four closed triangles are shown in Figure [31 The non-zero 
closure phases indicate the presence of an asymmetry in the disk. 



^http://kurucz. harvard.edu/grids. html 

^Cahbrated OIFITS data are available upon request. 
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3. Modeling the Disk of C Tau 

We fit a two-component geometric model to the MIRC data obtained for ( Tau. The 
model is compo sed of a uniform disk with an angular diameter of 0.40 mas {R=5.5 Rq] 
Gies et al.l 120071 ) to fit the central star and an elliptical, Gaussian surface brightness dis- 
tribution to model the circumstellar disk. To account for the asymmetry we detect in the 
closure phases, we modulated the ellipti cal Gaussian disk by a sinusoid as a func tion of 
projected azimuth or position angle (e.g. iMonnier et al.ll2006bt iThureau et al.ll2009l ). This 
creates a "skewed" disk model where the sinusoid causes the brightness distribution to peak 
on one side of the disk and places a depression in the brightness on the other side. The 
intensity distribution on the sky of the asymmetric, elliptical, Gaussian disk is given by the 
following functional form, 

y' = X sin 0maj + 1/ COS I 



X 



/disk = /o [1 + ^skew * COS*' 



XCOS( 



''skew 



'maj 



y sm I 



•'maj 



'maj 



exp < —4 In 2 



X 



+ 



^min / V ^maj 



(1) 

(2) 
(3) 



where 6'maj and 6. 



are the full-width at half maximum of the major and minor axes of the 
maj is the position angle of the major axis measured east of north, Askew is 



Gaussian disk, 

the amplitude of the sinusoidal modulation (0—1), 0skew is the position angle of the skew 
maximum intensity measured east of north, and Jq is the normalized central brightness of 
the disk. In this parametrization, +x runs in the direction of positive RA (east), +y in the 
direction of positive declination (north), and (p = arctan(x/?/). The shape of the asymmetry 
can be fiattened out to be more boxy in appearance by raising the sinusoid by the "skew 
power" p {p < 1.0). To avoid taking the root of a negative number, we take the absolute 
value of cos (0skcw — 0), raise it to the power p, and then multiply that value by the original 
sign of cos (0skew — 0)- Finally, we scale the star and disk contributions by the fraction of 
their if-band fiuxes (/star and /disk)- 

To fit the parameters of the ( Tau model we performed a Levenberg-Marquardt least- 
squares minimization using the IDL mpfilj^ rountine developed C. B. Markwardt. For each 
iteration, we computed an image of the geometric model and computed the visibilities and 
closure phases at the same u-v coordinates as the data and compared these with the observed 
values. Table [3] lists the model parameters derived for each epoch. We present the total 
reduced x^(all) for a solution and also break it down into separate contributions from the 
visibilities xli^"^) and closure phases x^(T3). The formal la uncertainties listed in Table [3] 



'http://cow.physics.wisc.edu/~craigm/idl/idl.htinl 
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were computed from the square root of the diagonal elements in the covariance matrix and 
scaled by xl- To test the reliability of the derived parameters, we fit the data from each 
night separately in addition to doing a global fit to data observed on multiple nights at the 
same epoch. For the 2007 Nov data, we decided to group only the 2007 Nov 14 (inner array) 
and 2007 Nov 19 (outer array) data in the global-epoch fit. We did this to ensure that the 
number of data points going into the global fit from the inner and outer arrays were roughly 
equal (to provide equal weights for each configuration). Additionally, the data from 2007 
Nov 14 yielded the best individual fit for the inner array; the data quality from 2007 Nov 11 
and 13 was not as good because of poorer seeing conditions. In Figures [2] and [3|, we overplot 
the model in each panel to compare with the measured visibilities and closure phases. Images 
of these models derived for 2007 Nov 14+19, 2008 Sep 26-28, 2008 Dec 10, and 2009 Nov 10 
are shown in Figure IH We note that the ( Tau data on 2008 Dec 10 were taken during poor 
seeing conditions that improved by the time we observed the second calibrator of the night. 
The changing seeing conditions will result in larger systematic calibration errors. We noted a 
trend of increasing system visibility on all baselines during that night, so we suspect that the 
calibration will affect the overall scaling of the disk size more than it affects the orientation 
or axis ratio of the disk. Additionally, the 2008 Dec 10 data only had one sampling in u-v 
space, where the other epochs and nights had several observations to improve confidence in 
the results. 

The largest source of systematic error lies in the calibration of the MIRC visibilities. The 
main contributions come from the photometric corrections for the amount of light in each 
fiber and errors in the angular diameters we assumed for the calibrator stars. We estimate 
that these effects can result in a ~10% systematic uncertainty in the squared visibilities. 
To investigate how this affects the model parameters, we generated 100 data sets for each 
epoch where the visibilities on each baseline were altered by a fraction of ±10%; the values 
of the offsets were drawn randomly from a Gaussian distribution. Because the visibilities are 
correlated across the 8 spectral channels, we varied all of the visibilities on a given baseline by 
the same fraction. Using these modified data sets, we derived a new set of best fit parameters 
for the disk model. We then estimated the systematic effects on a given parameter by taking 
the standard deviation of the parameter distribution generated from the modified data sets. 
TableHlshows the size of these effects on the fraction of flux in the disk, the major axis, minor 
axis, and position angle of the disk. These parameters are the most sensitive to the visibility 
calibration. The size of the systematic effects are typically 2-5 times the formal internal 
errors determined from the covariance matrix (Table [3]). The systematic errors tend to have 
a greater impact on the data from the inner array than the outer array. This is because a 
10% change in the visibility calibration will have a larger effect on shorter baselines where 
the source is less resolved and the visibilities are closer to 1.0. 
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Discussion 



4.1. Size and brightness of the disk in the H-hand 

Across the epochs with reliable visibility calibration (excluding 2008 Dec 10), the 
FWHM of the major axis of the disk ranges between 1.6-2.1 mas (0.20-0.27 AU). The range 
of variation does not significantly exceed the systematic uncertainties, so we cannot con- 
clusively state whether we measure a true change in the outer radius of the disk. The size 
in the i f-band is similar to the FWHM of the major axis of 1.8 mas measured in the K'- 
band by iGies et al.l (120071) an d is smaller than the 3.1-4.5 mas FWHM measured in Ha by 



Quirrenbach et al 



(Il997f ) and lTvcner et al.l (12004 ). 



It is interesting to note that the inner array data in 2007 Nov 11-14 tend to favor a 
puffier disk, where the ratio between the size of the minor axis relative to the major axis 
is larger than the ratio derived from the outer array data on 2007 Nov 19. If this effect is 
real, it might suggest a polar wind cona ponent along the direction of the minor axis (e.g. 
Kanaan et al.l l2008l : iKervella et al.l 120091 ) that is less resolved with the shorter baselines of 
the inner array and more resolved with the longer baselines of the outer array. We plan to 
follow-up on this in future observations; the newly installed MIRC photometric channels will 
also improve the visibility calibration and help to discern whether the effect is real. 

In general, a degeneracy exists between the size of a circumstellar disk and the ffux ratio 
of the star-to-disk contributions. This degeneracy breaks down if we successfully resolve the 
disk and can see the visibility curve ffatten out at the longest baselines, indicating that the 
disk is fully resolved and hence, the visibility amplitude is dominated by the contribution 
of the underlying unresolved star. As seen in Figure [2l the visibility amplitude begins to 
ffatten out on long baselines that sample projections along the major axis of the ( Tau 
disk, indicating that we are able to remove much of this degeneracy. From our MIRC 
observations, we ffnd that the star contributes on average about 0.55 ± 0.08 of the light in 



the i f- band. In comparison, the star contributes 0.41 of the ffux in the i^'-band (iGies et al. 



20071 ). iTouhami et al.l (120101 ) recently measured the near-IR excess ffux in ( Tau and other 
Be stars by assuming that the stellar ffux dominates in the visible part of the spectrum. The 
neglect of any disk ffux in the visible will result in an overestimate of the stellar ffux in the 
near-IR. Therefore, their upper limits on the ratio of stellar to total ffux of 0.76 in H and 
0.62 in K are consistent with the interferometric results. 
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4.2. Position Angle Variations and Disk Precession 



The position angle of the major axis of the ( Tau disk changes across the different epochs 
of MIRC observations. For comparison, we list the model parameters of prior results from 
the literature in Table [Sj there is some scatter among these position angle measurements 
also. We include here a measurement derived from CHARA Classic i^^'-band observations 
by Y. Touhami et al. (in preparation) that were made contemporaneously with the MIRC 
observations in 2008. These position angle variations probably reflect changes in the disk 
spatial flux distribution in the sky that are related to disk gas asymmetries. 



jKatolf 



T he disks of Be stars may develop a global, one-armed spiral 



[m 



2008 



1983 



Okazaki 



1991 



Oktariani &: Okazaki 



1 mode) instability 



9971: iPapaloizou et al.lll992l : iPapaloizou fc Savonijell2006l : lOgilvie 



20091 ). The oscillation mode appears as a spiral density wave that 



precesses prograde with the disk rotation and completes a 360° advance over a timescale 
of a few years. The progressive change in gas density with disk azimuth is probably the 
explanation for cyclic changes in the intensities of the v i olet V an d red R peaks of the 
Ha em ission profiles. Recent theoretical work by lOgilvid ( l2008l ) and lOktariani fc Okazaki 
( I2OO9I ) suggests that there is significant vertical motion associated with the spiral arm that 
might as s ume t he fo rm of a tilted disk that precesses with the one-armed oscillation mode. 
Hummell ( 1l998l ) and iHiratal (120071 ) both argued that precession of the disk tilt is the cause 
of long-term emission profile variations in some Be stars. 

The cyclic variations in the V/R ratio of the Ha emission peak strengths is well docu- 



mented in the case of C Tau ( iRivinius et al.ll2006l : iPoUmann &: Riviniusll2008l : iRuzdjak et al. 



2OO9I : IStefl et al.ll2009[ ). The evidence relating the V/R variations to the oscillation model 
is especially striking for the last three cycles when the Ha emission peaks varied with a 
cycle time of 1429 days and with the max imum occurring at a reference epoch of JD 
2,450,414 dstefl et all bood : ICarciofi et al]l2009h . The Ha variations follow a cycle where 
the violet side of the emission line peaks at oscillation phase t = (V > R), descends to 
V = R at T = 0.25 when a central, "shell", absorption component appears, then the red 
peak reaches a maximum at r = 0.5 {V < R), and as cends back thr o ugh V = R at t = 0.75 
where a third central emission peak feature appears. ICarciofi et al.l ( 120091 ) argue that these 
spectral variations are caused by the changing orientation of a spiral density enhancement 
around the disk azimuth. They also present evidence for an asymmetry in the disk based 
on the astrometric shift in the photocenter of the blue and red-shifted parts of the Br7 line 
measured with VLTI/Amber. 

We collected more recent Ha spectra to confirm that the V/R variation continued 
through the time span of the MIRC observations. The spectra were downloaded from the 
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BeSS database maintained at the GEPI laboratory of the Observatoire de Paris-MeudoiJ^ 
and they were augmented by two spectra from the University of Toledo Ritter Observatory 
(courtesy of E. Hesselbach and K. Bjorkman) and one from the Kitt Peak National Observa- 
tory Coude Feed telescope (courtesy of E. Grundstrom and V. McSwain). All these spectra 
have a resolving power of approximately R = 10, 000 or better. They were transformed 
to continuum flux normalized versions on a standard heliocentric wavelength grid. The in- 
tegrated Ha emission equivalent width (measured without correction for any photospheric 
component) was relatively constant over this t ime with a mean value of Wx = -15.2 A, sim- 
ilar to the mean over 1992 - 2008 of -15.5 A flStefl et al.ll2009h . We measured a simple V/R 
estimate as the ratio of the maximum flux excess above the continuum on the blue side com- 
pared to that on the red side of the Ha emission profile, and the results are shown in Figure [5] 
(together with the times of the MIRC observations). The V/R variation did indeed continue 
through this period, but the maximum occurred near JD 2,454,505 ±30 or app roximately 196 



Carciofi et al.l (120091) . This trend is 



Stefl et al.l (120091 ). For the purpose 



d earlier than predicted according to the ephemeris from 
consistent with the slow decrease in cycle time noted by 
of comparing the disk position angles from the MIRC results with the V/R phase r, we set 
r using the observed recent time of maximum but with the 1429 d cycle time unchanged. 
Note that there is evidence of a shorter, ~ 70 d varia tion in the measurements pr ior to the 
V/R maximum, first noted in earlier observations by lPoUmann &: Riviniud (120081 ) . that we 
will discuss below. A third, central peak sometimes appears in the Ha profile near phase 
r ~ 0.75 that might confuse the V/R estimate. The V/R measurements presented in Fig- 
ure 5 span the range from r = 0.91 to 0.48, and from a visual inspection of the Ha profiles, 
we only see evidence of a weak, central peak in six spectra obtained between HJD 2,454,433 
and 2,454,457 (r = 0.95 — 0.97). The stronger V and R peaks are well-separated from the 
low intensity, third peak in all six of these observations, so we doubt that the third peak has 
any significant influence on this particular set of V/R measurements. 

In Figure El we plot the position angle of the disk measured with MIRC and other 
published interfometric observations as a function of V/R phase r. We list the phase com- 
puted for each epoch in the bottom row of Table s [3] and [5l Note that we do not include 
in Figure E] the estimate for 1992 from the work of lOuirrenbach et al.l (119971 ) since the V/R 
variatio ns were not as well documented then and the cycle time was probably significantly 
longer (IRivinius et al.l l2006l : iRuzdjak et al.ll2009l ). It appears that the position angle did 
vary with r over the last cycle. A sinusoidal fit to the variation yields a mean position angle 
of < 0maj >— — 58?0 ± 1?4, a semiamplitude of 8?1 ± 1?7, and an epoch of maximum position 
angle at r = 0.23 ± 0.03. The residuals of the fit are reduced from a rms = 4?6 for a simple 



^http://basebe. obspm.fr/basebe/ 
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mean to a rms = 2?2 for the sinusoidal fit, and an F-test indicates that such a reduction 
would only occur 3% of the time for random errors. Thus, we suggest that the position 
angle variations are significant and are probably related to the V j R variation. We note that 
the fitted mean position angle from inte rferometry agrees well with the mean from linear 
polarization measurements, -58? 1 ± 1?2 jMcDavidI Eoogl : Istefl et al.lboOQh . 



We suspect that the position angle variations result from a tilt of the disk that is 
associated with the one-armed spiral density enhancement. A cartoon model for the variation 
is presented in Figure [7] that shows the orientation of the circumstellar disk and the density 
enhancement over the V/R cycle. For the sake of simplicity, the tilt is shown as constant 
at all disk radii (for a planar geometry), although we suspect that in reality the disk tilt 
is largely confined to radii where the density enhancement is largest. The geometry of the 
projected disk in the sky is described by four parameters: the mean position angle of the 
disk normal ao measured east from north, the average disk inclination < i > (equal to the 
inclination of the spin axis of the Be star), the precession cone semiangle 6, and the time 
variable azimuth of the precession axis iP{t) (measured relative to the line of sight). In the 
thin disk approximation, the projected position angle of the disk major axis will vary as 



V) = ao + 90° — arctan(sin'?/'(r) tan^/ sin < i >] 



and the projected ratio of the minor to major axes will be related to the time variable 
inclination by 

cos i{r) = cosiP{t) sin 6' sin < i > + cos 9 cos < i > . 

We assume that the disk tilt follows the prograde precession of the one-armed spiral mode 
and th at the disk approaches us in the no rth-west sector (based upon the interferometric 



results; Vakili et al. 



1998 



Stefl et al.l 120091 ). Then, we derive from the sinusoidal fit of the 
148?0 ± 1?4, 9 = 8°l± 1?7, and a temporal relation for the 
= 360°(r+ (0.52±0.03)). Finally we set < i >= 92?8 based 
0.57) = 85° (angle between the disk normal directed 



position angle variation, = - 
azimuthal precession angle iP{t) 
upon the apparent inclination of i(r 
south-west and the line of sight) derived by ICarciofi et al.l (120091 ) from their disk model of 
the VLTI/ Amber observations. 

We next consider the relationship between the disk tilt geometry and the location of 
the density enhancement. We place the enhancement maximum in the tilted disk plane at 
an arbitrary distance from the Be star of 2.5Ri, (Carciofi et al. 2009), and it is marked as 
a gray circle in each panel of Figure 7. Recall that r = is defined by the time of V/R 
maximum when the density enhancement is located in the approaching part of the disk in 
the plane of the sky. This relation defines the azimuthal placement of the enhancement. In 
the cartoon model, the precession motion brings this enhancement almost directly in front 
of the Be star at r = 0.25, which is presumably the reason for the strong, shell absorption 
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feature that appears in the Ha profile then. Half a cycle later at r = 0.75 the density 
enhancement is occulted by the Be star, but we speculate that the spiral arm extension is 
large enough that some high density regions are still visible and contribute to the appearance 
of a third emission peak near the center of the Ha profile. The azimuthal relation we find 
for '?/'(r) places the density enhancement maximum near the line of the nodes between the 
tilted disk plane and equatorial plane of the Be star, and it will be interesting to see if such 
a placement will be found in future three dimensional models for the one-armed oscillation. 
Note that for this location and disk inclination angle, we expect that the photocenter of the 
enhancement will always be found along the line of the projected stellar equator (dotted line 
in Fig. 7). Thus, the observed position angle variations result not from the shifting position 
of the enhancement but from the precession of the extended disk whose tilt is presumably 
generated by vertical motions associated with the enhancement. Our simple model also 
predicts that the observed ratio of disk minor to major axes will be smallest around r = 0.28 
and 0.68, which is consistent with our finding of the smallest ratio in the 2008 Dec 10 data 
at V/R phase r = 0.21. 



The fact that C, Tau has a binary companion in a 133 d orbit (IRuzdjak et al.l 120091 ) 
probably means that any tilt of the Be star's disk will be modulated by the tidal force of 
the binary companion. Approximately twice each orbit, a tilted disk will experience a tidal 
torque in the direction of coalignment with the orbital plane, and this results in a nodding 
motion that is seen, for exarn ple, in the precessing disk and jets of massive X-ray binary 
SS 433 (ICoUins fc Scherll2002[ ). For prograde precession, the nodding period is 



Pn 



1 PpPh 

2P„-P 



where Pp and Pb are the precessional and binary periods, respectively. The predicted nodding 
period, Pn = 73.1 d, is quite close t o the observed V/R modulation period of 69.3 ± 0.2 d 
discovered by lPoUmann &: Riviniud ( 120081 ) and which is evident in the recent V/R variations 

(Fig. ED. 

A lingering difficulty for the precessing disk model is the fact that the intrinsic polariza- 
tion angle has remained remarkably stable over the last decade ( lMcDavidlll999l : IStefl et al. 
20091 ). showing no evidence of changes as large as those seen in the interferometric data. We 
suspect that this difference arises in the radial dependence of the disk shape. The polariza- 
tion measurements probe scattering radiation from the innermost part of the disk where the 
disk is probably coaligned with the stellar equator (the probable source of entering disk gas) . 
The tilt, on the other hand, is associated with vertical motions produced by the one-armed 
spiral oscillation that attain maximum amplitude at several stellar radii out into the disk. 
Thus, we suggest that the sensitivity of the interferometric observations to disk emission 
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at larger radii is the reason why the tih oscillation is detected by interferometry and not 
detected by polarization measurements that reflect conditions in the inner disk. 

4.3. Nature of the Asymmetry in the Light Distribution 

The non-zero closure phases we measure in ( Tau imply the presence of an asymmetry 
in the disk and/or the star light distribution. In general, the position angle of the peak 
of the skewed distribution lies within ~10° of being perpendicular to the major axis of the 
disk. However, because there is a fairly sharp transition between the "bright" and "dark" 
sides of the disk, a small change in the angle of the asymmetry has a large effect on how the 
transition line intersects the outer edges of the disk and changes the brightness distribution 
there. 

From the model images in Figure HJ it appears that in 2007 Nov, the upper (north- 
eastern) half of the disk is brighter than the lower (south-western) half. In 2008 Sep, the 
north-western side of the major axis appears brighter than the south-eastern side. In 2009 
Nov, the south-eastern side of the major axis is brighter than the north-western side. These 
locations agree roughly with the quadrant where we expect to see the density enhancement 
based on the precession mode (gray circle in Fig. [Tj), provided that the spiral arm extension 
covers a larger area of the disk than marked in the cartoon model. A more sophisticated test 
of the one-armed spiral oscillation model would require adopting a more realistic brightness 
distribution for the density enhancement in addition to taking into account how the Be star 
is occulted by the disk. 

In 2008 Dec, the disk appears very thin. The skew parameters (Askew = 1-0, p = 0.001) 
indicate that the upper side of this thin disk is totally bright and the lower side is totally 
dark. This suggests that the disk is tilted slightly toward the upper half of the star, similar 
to a phase between r=0.25 and r=0.375 displayed in the panels of Figure [71 Given the 
uncertainties in measuring the cycle length and the time of V/ R maximum, this is roughly 
consistent with our estimate of r=0.21 during this epoch. The model image for 2008 Dec in 
Figure m shows a slight enhancement in brightness toward the north-west side of the major 
axis. However, we note that the skew axis 0skcw = 44?11 ± 0?09 deg is almost identical to 
the disk normal at a position angle 90° + 0disk = 44?40 ± 0?10 deg, indicating that we cannot 
conclusively determine whether the brightness is skewed toward one side of the major axis. 

If the binary companion of ( Tau were bright enough to effect our observations, we 
would expect to see a periodic variation in the closure phases; the amplitude of the variation 
revealing the flux ratio and the periodicity corresponding to the separation on the sky (e.g. 
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Monnierl 120071 ) ■ In Table O we list the predicted location of the binary companion at the 
times of the MIRC observations based on estimates of the orbital parameters. To compute 
th e separation (ph ^ ■n) and position angle (^bin), we used the spectroscopic parameters derived 



bv lRuzdiak et al.l (120091 ) in their Table 2, Solution 2 (P = 132.987 days, = 7 A3 kms-\ e 



= 0, a; = 0). Assuming that the binary orbit lies in the same plane as the disk and that the 
epoch of maximum radial velocity (TpiVmax = 2,447,025.6 HJD) occurs when the companion 
is approaching from the north-west, we adopted an orbital inclination of i = 92?8 and a 
position angle of the li ne of nodes of fl = — 58?0. If we assume that the B2 Ilipe primary 
has a mass of 11.2 Mq (IGies et al.l 120071 ). these orbital parameters yield a secondary with a 
mass of 0.94 Mq. Using Kepler's Third Law, these mass estimates imply a semi- major axis 
for the binary orbit of 1.17 AU, or a = 9.17 mas at a distance of 128 pc. If the companion 
is a main sequence star, the secondary mass corresponds roughly to a G4 spectral ty pe with 
an absolute magnitude of My = 5.0, V - K = 1.5, and K - H = 0.1 dCro] boool ) . Using 
the distance modulus implied by the Hipparcos parallax, this corresponds to an apparent 
if-band magnitude of 9.0. Compared with the 2MASS magnitude of ( Tau (if=3.0), this 
corresponds to a magnitude difference of AH = 5.9. 

It could be argued that we measure a small periodic amplitude variation in the closure 
phases plotted in Figure |3l However, the data quality is not sufficient to compute a full 
binary fit. Adding a binary component to our skewed elliptical Gaussian models at a fixed 
fiux ratio and a position given by the values in Table Ej does not improve the significantly. 
Even though we do not detect the companion, the variation in the closure phase residuals 
allows us to put an upper limit on the binary fiux ratio. The maximum deviation from 
the best fit skewed disk model ranges from 2?2 in 2007 Nov to 7?1 in 2009 Dec. If this 
variation is due only to the signature of a binary companion, then we can place a lower 
limit on the if-band magnitude difference between ( Tau and the companion of 5.3 mag and 
3.3 mag, respectively. These values are consistent with the expected magnitude difference 
of AH = 5.9 estimated from the orbital parameters we assumed. The observed magnitude 
limit has interesting implications for the possibility of a hot subdwarf companion. If we take 
the effecti ve temperature an d radius for the hot subdwarfs in Per (IGies et al.lll998[ ) and 
FY CMa ( iPeters et al.l 120081 ) and assume a Planck flux distribution, then we would predict 
magnitude differences of AH = 2.5 (0 Per subdwarf) and 4.4 (FY CMa subdwarf), which 
would probably be detectable. So if the companion is a subdwarf, it is probably cooler 
and/or smaller than the subdwarfs detected in the other two cases. 
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Conclusions 



We obtained four epochs of interferometric measurements on tlie Be star ( Tau using 
the MIRC beam combiner at the CHARA Array. By fitting the disk with a skewed elliptical 
Gaussian model, we determine a full width at half maximum size of the major axis of ~ 
1.8 mas and estimate that the central Be star contributes 55% of the light in the H-hand. 
Combining our results with previous interferometric measurements, we observe a change in 
the position angle of the disk over time. A correlation between the position angle and the 
V/ R phase of the Ha emission line variation suggests that the tilt of the disk around ( Tau 
is precessing. The tilt could be gerierated by vertical motions o f the gas caused by the spiral 
density enhancement (jOgilviell2008t lOktariani fc Okazakill2009l ) as it moves through the disk. 
We also measure an asymmetry in the light distribution of the disk that roughly corresponds 
to the expected location of the density enhancement in the spiral oscillation model. 

We plan to continue monitoring changes in the structure and orientation of the disk of 
( Tau with future observations at the CHARA Array. Ultimately, these observations can 
be used to test predictions from one-armed spiral oscillation models (e.g. Berio et al. 1999; 
Meilland et al. 2007; Jones et al. 2008; Carciofi et al. 2009). 
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Table 1. CHARA-MIRC Observing Log for C Tau 



UT Date Configuration Calibrators 



2007 


Nov 11 


S2- 


-E2- 


-Wl- 


-W2 


C Per 


2007 


Nov 13 


S2- 


-E2- 


-Wl- 


-W2 


C Per 


2007 


Nov 14 


S2- 


-E2- 


-Wl- 


-W2 


0- Cyg, C Per 


2007 


Nov 19 


Sl- 


-El- 


-Wl- 


-W2 


C Per, 10 Aur 


2008 


Sep 26 


Sl- 


-El- 


-Wl- 


-W2 


( Per, 6 Gem 


2008 


Sep 27 


Sl- 


■El- 


-Wl- 


-W2 


C Per, 9 Gem 


2008 


Sep 28 


Sl- 


■El- 


-Wl- 


-W2 


CPer 


2008 


Dec 10 


Sl- 


■El- 


-Wl- 


-W2 


C Per, 15 LMi 


2009 


Nov 10 


81- 


-El- 


-\Y1- 


-\Y2 


HR 185. C Pel-. 71 Ori. 79 Cnc 



Table 2. Calibrator Diameters 



Calibrator 


HD Number 


SpT 


V 


H 




Ref 


E(B - V) 


R 


LD Diameter 








(mag) 


(mag) 


(K) 




(mag) 




(mas) 


HR 485 


HD 10348 


KO HI 


5.97 


3.83 


4885 


1 


0.046 ± 0.031 


3.1 


0.879 ± 0.046 


C Per 


HD 24398 


Bl lab 


2.87 


2.62 


21950 


2 


0.355 ± 0.012 


2.88 ± 0.09 


0.645 ± 0.026 


10 Aur 


HD 32630 


B3 V 


3.16 


3.76 


16600 


3 


0.024 ± 0.007 


3.1 


0.444 ± 0.012 


71 Ori 


HD 43042 


F6 V 


5.20 


3.83 


6485 


4 


± 0.007 


3.1 


0.597 ± 0.021 


9 Gem 


HD 50019 


A3 III 


3.60 


3.23 


8300 


5 


0.033 ± 0.009 


3.1 


0.796 ± 0.022 


79 Cue 


HD 78715 


G5 III 


6.0 


4.09 


5050 


6 


± 0.007 


3.1 


0.747 ± 0.062 


15 LMi 


HD 84737 


GO V 


5.08 


3.61 


5830 


7 


0.046 ± 0.024 


3.1 


0.849 ± 0.022 


<y Cyg 


HD 202850 


B9 lab 


4.25 


3.86 


11000 


8 


0.237 ± 0.009 


3.1 


0.574 ± 0.017 



References. — (1) Hekker & Melendez 2007; (2) Huang & Gies 2008; (3) Lyubimkov et al. 2002; (4) Lambert & Reddy 2004; 
(5) Malagnini et aL 1982; (6) de Laverny et aL 2003; (7) Ramirez et al. 2007; (8) Markova & Puis 2008 



Table 3. Disk Modeled as an Elliptical Gaussian Modulated by a Sinusoid 



Parameter 


2007Novll,13 


2007Novl4 


2007Novl9 


2007Novl4,19 


2008Sep26-28 


2008DeclO 


2009NovlO 




S2E2W1W2 


S2E2W1W2 


S1E1W1W2 


S2E2W1W2 
S1E1W1W2 


S1E1W1W2 


S1E1W1W2 


S1E1W1W2 


/disk 


0.5948±0.0073 


0.5464±0.0078 


0.5036±0.0065 


0.5114±0.0065 


0.4717±0.0060 


0.4669±0.0243 


0.3480±0.0051 


/star 


0.4052±0.0073 


0.4536±0.0078 


0.4964±0.0065 


0.4886±0.0065 


0.5283±0.0060 


0.5331±0.0243 


0.6520±0.0051 


6'inaj (mas) 


1.662±0.042 


1.527±0.045 


1.410±0.025 


1.599±0.036 


1.633±0.031 


1.116±0.061 


2.057±0.060 


6'min (mas) 


0.691±0.027 


0.598±0.047 


0.308±0.015 


0.331±0.023 


0.401±0.010 


0.020 


0.346±0.016 




0.416±0.019 


0.392±0.033 


0.218±0.011 


0.207±0.015 


0.246±0.008 


0.018±0.001 


0.168±0.009 


0disk(°) 


-61.82±0.86 


-62.07±0.75 


-58.83±0.71 


-60.92±0.66 


-51.55±0.51 


-45.60±0.10 


-57.33±0.86 


-^skew 


0.92±0.11 


0.35±0.70 


0.34±0.02 


0.40±0.08 


0.49±0.09 


1.00 


0.21±0.02 


(pskew{°) 


29.49± 0.65 


29.03± 0.22 


32.42± 0.37 


30.44± 0.62 


32.88± 0.53 


44.11± 0.09 


41.28± 1.11 


P 


1.00 


0.05±1.35 


0.001 


0.16±0.09 


0.52±0.11 


0.001 


0.001 


X^(all) 
xl (T3) 


2.25 


1.11 


1.12 


2.16 


1.27 


1.16 


1.33 


3.40 


1.31 


1.09 


2.46 


1.09 


1.74 


1.47 


0.64 


0.94 


1.33 


1.84 


1.61 


0.58 


1.16 


HJD - 2,400,000 


54417.0 


54419.0 


54423.8 


54421.4 


54736.9 


54810.8 


55145.9 


T {V/R Phase) 








0.942 


0.162 


0.214 


0.449 



Note. — 



Values without error bars indicate that the fitting routine reached the limit of the search range. 
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Table 4. Systematic Errors in Model Parameters Estimated Through a Monte Carlo 

Analysis 



Parameter 


2007Novll-14 


2007Novl4,19 


2008Sep26-28 


2008DeclO 


2009NovlO 




S2E2W1W2 


S1E1W1W2 


S1E1W1W2 


S1E1W1W2 


S1E1W1W2 






S2E2W1W2 








^disk 


±0.027 


±0.020 


±0.030 


±0.091 


±0.027 


0maj (mas) 


±0.40 


±0.12 


±0.18 


±0.15 


±0.45 


^min (mas) 


±0.15 


±0.074 


±0.047 


±0.080 


±0.073 


0,U,-,k(") 


±1. •-!.() 


±2.4 


±2.<S 


±2..") 


±4.:-! 


Table 5. Previous Interferometric Results 


Parameter 


1992.82 


1999.16 


2005.93 


2006.95 


2008.80 


/star 


0.30 


0.814±0.012 


0.414±0.029 






6'maj (mas) 


4.53±0.52 


3.14±0.21 


1.79±0.07 






^min/ ^maj 


0.28±0.28 


0.310±0.072 


0.09±0.22 






0disk (°) 


-58±4 


-62.3±4.4 


-52.2±1.7 


-58±5 


-51.8±4.0 


Observatory 


Mklll 


NPOI CHARA-Classic 


VLTI-AMBER 


CHARA-Classic 


Filter 


Ha 


Ha 


K' 


K 


K' 


HJD - 2,400,000 48915.8 


51238.8 


53709.3 


54081.7 


54758.9 


r {V/R Phase) 




0.577 


0.306 


0.567 


0.178 


Reference 


1 


2 


3 


4,5 


6 



References. — (1) Quirrenbach et al. 1997; (2) Tycner et al. 2004; (3) Gies et al. 2007; (4) Stefl et 
al. 2009; (5) Carciofi et al. 2009; (6) Touhami et al., in preparation 



Table 6. Predicted Location of the Binary Companion During the MIRC Observations 



2007Novl4 2007Novl9 2008Sep26-28 2008DeclO 2009NovlO 



HJD - 2,400,000 54419.0 54423.8 54736.9 54810.8 55145.9 

Pbin (mas) 7.59 6.24 9.13 8.86 8.51 

<pun (°) 120.1 119.0 302.3 121.3 300.9 
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200 
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-100 



-200 
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II 



2007NOV11-I-13 



200 100 -100 -200 
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200 
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-100 



-200 




200 100 -100 -200 
u (10* cycles) 



200 



^ 100 



■100 



-200 




200 100 -100 -200 
u (10* cycles) 



200 



100 



■100 



-200 



■I I ■ 



\ 



-I- 



\ 



2008Decl0 



200 100 -100 -200 
u (10* cycles) 



200 



100 



10 

u 
u 



■100 



-200 




200 100 -100 -200 
u (10* cycles) 



Fig. 1. — u-v coverage on the sky during the MIRC observations of C Tau in 2007-2009. 
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100 150 200 
B/X (xlO') 



50 100 150 200 
B/A (xlO') 



Fig. 2. — Squared visibilities measured for C Tau using MIRC in 2007-2009. For clarity, the 
measurements during each epoch are grouped by proximity in position angle of the baseline. 
The small inset panels show the projection of observed u-v points on the plane of the sky. 
The solid lines show the best global-epoch fits given in Table [3l There are a number of lines 
in each plot to show the model at each of the observed u-v projections. 
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S2-W1 -W2 


S2-E2-W1 






S2-E2-W2 

Ai J_ii_* JA 


E2-W1 -W2 


2007Nov14 


7»T? f » ? ? ' 


50 1.55 1.60 1.65 1.70 1.75 1.55 1.60 1.65 1.70 1.75 1. 
Wovelength (/xm) Wovelength (/xm) 


S1-E1-W1 


S1-E1-W2 






S1 -Wl -W2 
2008Sep26-28 


El -Wl -W2 





S1-E1-W1 


S1-E1-W2 




^ ^ i, ± * H 






S1 -Wl -W2 
' i * * i I A A 


El -Wl -W2 


2007Nov19 






50 1.55 1.60 1.65 1.70 1.75 1.55 1.60 1.65 1.70 1.75 1. 
Wovelength (/zm) Wovelength (/zm) 


S1-E1-W1 


S1-E1-W2 










S1 -Wl -W2 


El -Wl -W2 

t-UitiJ-i-: 


2008DeclO 
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Fig. 3. — Closure phases measured on the four closed triangles during our MIRC observations 
in 2007-2009. The solid lines show the best global-epoch fits given in Table [3l The multiple 
lines show how the model changes across the observed u-v projections. 
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Fig. 4. — Best-fit geometric models determined for ^ Tau during the epochs of the MIRC 
observations. The spectroscopic V/ R phase r is indicated in the bottom right of each paneL 
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Fig. 5. — The logarithm (base 10) of the V/R ratio of the Ho; emission peaks as a function 
of time (hehocentric Juhan date). The vertical tick marks at bottom indicate the times of 
the MIRC observations. 
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Fig. 6. — Disk long axis position angle as derived from interferometry plotted against V/R 
phase T. The measurements are repeated over two cycles to emphasize phase continuity. 
The filled black squares represent our MIRC observations (Table 3), the gray squares are the 
previously published near-IR measurements listed Table 5, and the open diamond indicates 
the Ha result from Tycner et al. (2004). The solid line shows a sinusoidal weighted fit of the 
variation, and the dotted line indicates the mean position angle determined independently 
from linear polarization observations. 
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Fig. 7. — Cartoon depiction of the disk precession variations as seen in the sky. The Be star 
is shown as a black circle, the circumstellar disk as a gray ellipsoid, and the mean position 
angle of the long axis of the disk (aligned with the stellar equator) as dotted lines extending 
from the disk. The panels show the change in the precession of the disk tilt and of the 
one-armed spiral density maximum {small gray dot) for eight steps in the precession cycle 
from T — 0, the time of V/R maximum. The disk rotation and precession both advance 
around the Be star's spin axis (pointed to the south-west) so that the north-western part of 
the disk approaches us. 



